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u. ABSTRACT 

(~| ■ We present an observational study about the impacts of the interactions in the kinematics, stel- 

PLhI lar populations, and oxygen abundances of the components of the galaxy pair AM 2322-821. 

A fairly symmetric rotation curve for the companion (AM 2322 B) galaxy with a deprojected 
velocity amplitude of 110 km s"' was obtained, and a dynamical mass of 1 . 1 - 1 .3 x 1O'°M0 
within a radius of 4 kpc was estimated using this deprojected velocity. Asymmetries in the 
^ ■ radial velocity field were detected for the companion, very likely due the interaction between 

the galaxies. The interaction between the main and companion galaxies was modelled us- 
ing numerical N-body/hydrodynamical simulations, with the result indicating that the cur- 
rent stage of the system would be about 90 Myr after perigalacticum. The spatial varia- 
tion in the distribution of the stellar-population components in both galaxies was analysed 
using the stellar population synthesis code STARLIGHT. The companion galaxy is domi- 
{Jr^ ' nated by a very young (t < 1 x 10^ yr) population, with the fraction of this population to 

, the total flux at /1 5 870 A, increasing outwards in the galaxy disc. On the other hand, the 

stellar population of AM 2322A is heterogeneous along the slit positions observed. Spatial 
^-s I profiles of oxygen abundance in the gaseous phase were obtained using two diagnostic dia- 

' grams(/?23=([Oii]^3727-i-[Ora]^4959-H[Om]^5007)/H/3vs. [Oiii]^5007/[Oii]/l 3727 and 

[O iii]^ 5007/[O ii]^ 3727 vs. [N ii]^ 6584/[0 ii]^ 3727), where we compared the observed val- 
ues with the ones obtained from photoionization models. Such gradients of oxygen abundance 
are significantly flatter for this pair of galaxies than in typical isolated spiral galaxies. This 
k> ' metallicity distribution is interpreted as the gradients having been destroyed by interaction- 

^ , induced gas flows from the outer parts to the centre of the galaxy. 
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1 INTRODUCTION 

Galaxy interactions and merger events play an important role 
on the evolution and the stellar formation history of galax- 
ies. Interacting/merging galaxies show enhanced star forma- 
tion when com pared with isola ted obj ects, as indicated by dif- 
ferent studies (Kennicutt et al. 1987; Sekiguchi & Wolstencroft' 



19921: IPonz elli & Pastoriza 1997: Barton Gillespie et al. 2003; 
Geller et al i booel; IWoods & Gellerll2007D . This enhancement has 
been observ ed as being a function of the projected galaxy pair sep- 
aration (e.g. iBarton et al.ll2000l: iLambas et aDl2003l; InSoUc et al.1 



l2004h . as we ll as being stronger in low-mass than in high-mass 
galaxies (e.g. IWoods & GelleJl2007l ; lEriison et al.l20(j^ . 
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The induced star formation associated with the gas motions 
created by the interaction also is expected to have an impact in the 
chemical state of the galaxies. As star formation enriches the inter- 
stellar medium via nucleosynthesis, inflows of metal poor gas from 
the outer parts of the galaxy can decrease the metallicity in inner re- 
gions and modify the radial abundance gradients across spiral discs. 
In fact, studies have found that interacting galaxies do not follow 
the well established correlation between luminosity and metallicity 
found in normal disc galaxies. The central regions of these galax- 
ies are u nderabundant when compared to isolate d galaxies of simi- 
lar mass jKewlev et al. 2006; Rupke et al. 2008; Ellison et alj2008l; 
iMichel-Dansac et alJl2008l ; IPeeples et alil2009l) . 
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Shallower metallicity gradients have been found in barred 
galaxies and explained by the action of inward and outward radial 
flows of interste llar gas i nduced by the non-axisymmetric poten- 



tial of bars (e g. iFriedli et al. 1994; Sellwood & Wilkinson 1993; 



iRov&Walsd f 19971) . N-body/Smoothed Particle Hydrodynamics 
numerical simulations of equal-mass mergers, although ignore the 
presence of ongoing star formation, predict that the radial metallic- 
ity gradients in disc galaxies flatten shortly after th e first pericen - 
tre passa ge, due to the radial mixing of gas ( RupkeetaLll2010 l). 
Recently, iKewlev et all j201Cl) deteiTnined the metallicity gradi- 
ents for 8 galaxy pairs with mass ratio near unity and showed 
that they were signifi cantly shallower than those in isolated spiral 
galaxies. Futhermore, iKrabbe et al.l ( 1200 Sh found for the interact- 
ing pair AM 2306-721 with a mass ratio of 2: 1 that the disc of the 
main galaxy showed a clear radial metallicity gradient, while the 
secondary presented a relatively homogeneous oxygen abundance. 
These authors interpreted the absence of abundance gradient in the 
secondary galaxy in terms of mixing the low metallicity gas from 
the outer parts with the rich metal gas of the centre of the galaxy. Is 
this picture always reproduced in others minor mergers? This ques- 
tion was not adressed before by theoretical merger simulations and 
observations, except the above study of AM 2306-721. Therefore, 
more observational studies of mergers with different mass ratios 
should provide use ful insights to answer the a bove question. We 
have selected from lperreiro & Pastorizal (|2004 several systems to 
study the effects of the kinematics, stellar population and gradient 
abundances of the galaxies in minor mergers, where t he first results 
of thi s programme were presented for AM 2306-721 dKrabbe et al.l 
l2008h . 

This paper presents the results for the system AM 2322-821, 
which is morphologically very similar to AM 2306-721, but its 
mass ratio is much lower than the latter. AM 2322-821 is composed 
of a SA(r)c galaxy with disturbed arms (hereafter, AM 2322A) in 
interaction with an irregular galaxy (hereafter, AM2322B). Both 
galaxies contain very luminous H ii regions with Ho- luminosity in 
the range of 2.53 X lO'' <L(Hq-) < 1.45 x lO^'erg s"' as esti- 
mated from Hff ima ges and high star for mation rate in the range of 
0.02 to 1.15 Mo/yr frerreiro et al.lbOOSk . 

The present paper is organized as follows: in Section |2] we 
summarize the observations and data reduction. The gas kinemat- 
ics of each galaxy and the numerical N-body/hydrodynamical sim- 
ulations of the interaction are presented in Sections|3]and|4] respec- 
tively. In Section[5] we present the stellar population synthesis. The 
metallicity analysis is in Section[6l and the conclusions are summa- 
rized in Section|7] 



2 OBSERVATIONS AND DATA REDUCTION 

Long slit spectroscopic data were obtained on 29/30 June 2006, 
01/02 July 2006, and 27/28 July 2008 with the Gemini Multi- 
Object Spectrograph at Gemini South, as part of poor weather pro- 
grammes GS-2006A-DD-6 and GS-2008A-Q-206. Spectra in the 
range 3 450 to 7 130 A were acquired with two settings with the 
B600 grating, and the 1" slit, keeping a compromise between spec- 
tral resolution (~ 5.5 A), spectral coverage and slit losses (due to 
the Image Quality = ANY constraint). The blue setting provided 
a wavelength coverage of 3 450 to 6 280A and the red setting of 
4 280 to 7 130A at about the same spectral resolution. The frames 
were binned on-chip by 4 and 2 pixels in the spatial and spectra 
directions, respectively, resulting in a spatial scale of 0.288 " pxr' , 
and0.9Apxr' dispersion. 



Spectra were taken at four different position angles on the sky, 
with the goal of observing the nucleus and the brigthest regions 
of the galaxies, as well as, one spiral arm that is away from the 
main galaxy. PA=59° is in the slit position crossing the nucleus of 
AM 2322A; PA=28° is cutting across main body of primary, but not 
across the nucleus (offset of about 8" NW from nucleus); the slit 
position at PA=60° is located off the disc of the main component, 
along the NW spiral arm (located between the main and secondary 
component). The PA=318° slit position is oriented along the main 
axis disc of the secondary component (AM 2322B) and also along 
the AM 2322A NE spiral arm. The PA=60° was observed only in 
the red spectra, and thus in this slit, it was possible only to study 
the ionized gas kinematics, but not their stellar population and the 
gas phase 0/H abundance. Figure[T]shows the four slit positions on 
the GMOS-S r' acquisition image. 

The exposure time on each single frame was limited to 700 
seconds to minimize the effects of cosmic rays, with multiple 
frames being obtained for each slit position to achieve a suitable 
signal. The slit positions are shown in Fig.[I] superimposed on the 
r'-band image of the pair. Table[T]gives the journal of observations. 
Conditions during both runs were not photometric, with thin cirrus 
and image quality in the range 0.5 " to 2.0 " (as measured from 
stars in the acquisition images taken just prior to the spectroscopic 
observations). 

The spectroscopic data reduction was carried out using the 
GEMINI. GMOS package as well as generic IRAF tasks. We followed 
the standard procedure: (1) the data were bias subtracted and flat- 
fielded; (2) the wavelengh calibration was established from the Cu- 
Ar arc frames with typical residuals of 0.2 A and applied to the 
object frames; (3) the individual spectra of same slit positions and 
wavelengh range were averaged with cosmic ray rejection; (4) the 
object frames were sky subtracted interactively using the gsskysub 
task, which use a background sample of off-object areas to fit a 
function to the specified rows, and this fit is then subtracted from 
the column of each spectra; (5) the spectra were relative flux cali- 
brated using observations of a flux standard star taken with the same 
set up as the science observations; (6) finally, one-dimensional 
spectra were extracted from the two-dimensional spectra by sum- 
ming over six rows along the spatial direction. Each spectrum there- 
fore comprises the flux contained in an aperture of 1" x 1.73". As- 
suming a distance of 49.6 Mpc for AM 2322-821 system barycenter 
(estimated from the radial velocities derived in Section |3]and total 
masses of AM2322A and AM2322B, estimated in SectionO and 
adopting /fo=75kms"' Mpc"'), this aperture corresponds to a re- 
gion of 241x419 pc- for AM 2322-821. The nominal centre of each 
slit was chosen to be the continuum peak at A 5735A. The mismatch 
between the blue and red spectra for all the apertures extracted was 
lower than about 5 % over the spectral range where both spectra 
overlapped. We decided not to co-add them. 



3 IONIZED GAS KINEMATICS 

The radial velocity was estimated from the strongest emission lines 
present in the spectra, namely i4861,[Oin] ^5007, Ha /16563 
and [Nn] /16584. The final radial velocity for each spectrum was 
obtained by averaging the individual measurements from the de- 
tected emission lines, and the errors were estimated from the stan- 
dard deviation of the individual measurements around the mean. 

The inclination of each galaxy with respect to the plane of the 
sky was computed as cos(0 = bja, where a and b is the minor and 
major semi-axes of the galaxy, respectively. The minor and major 
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Figure 1. GMOS-S r'-band image of AM 2322-821 with tlie observed slit 
positions. 



Table 1. Journal of observations 
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i I f 



Date (UT) 


Exposure time(s) 


PAO 


Ai(A) 


2006/06/29 


3x600 


318 


4280-7130 


2006/06/30 


6x600 


60 


4280-7130 


2006/07/01 


2x600 


28 


4280-7130 


2006/07/01 


3x700 


59 


4280-7130 


2006/07/01 


3x600 


318 


4280-7130 


2008/07/27 


6x600 


28 


3450-6270 


2008/07/27 


6x600 


59 


3450-6270 


2008/07/27 


6x600 


318 


3450-6270 



Figure 3. Same asg] but for AM 2322B. 



semi-axes as well as, the position angle of the major axis of each 
galaxy were obtained from the acquisition images in the r' filter, 
using a simple isophotal fitting with the IRAF stsdas.ellipse task. 
The fitting results for the position angle of the major axis and the in- 
clination as a function of the projected distance in arcseconds along 
the isophotal major axis are shown in Figs|2]and|3] for AM 2322A 
and AM2322B, respectively. 

For AM2322B, the position angle of the major axis and the 
inclination of the galaxy are nearly constant out to a 5" radius. 
The resulting values for the disc inclination and position angle of 
the line of nodes are i = 63° and tf/o = 318°, respectively. For 
AM 2322 A, the position angle of the major axis and the inclination 
of the galaxy as measured from the isophotal fitting show signif- 




Figure 2. Results of the isophote fitting as a function of the distance along 
the major axis a in arcseconds for AM 2322A. Top panel: inclination of the 
galaxy. Bottom panel: position angle of the major axis. 



leant variations with radius. According to iBender & Moellenhofil 
( Il987l) , the interaction among galaxies can cause isophotal twist- 
ing, so these variations can be due to the disturbed morphology of 
the galaxy, added to the fact that the GMOS r'-band image is not 
deep and star forming regions and spiral structure are dominating 
in these wavelengths. If we considered the external isophotes of 
AM 2322A, the values found are about tf/Q = 75° and ; = 20°. 

iFerreiro & Pastorizal (|2004|) estimated the inclination of the 
galaxies in this pair measuring the major and minor diameters of the 
24 mag arcsec"^ isophote from images in the B filter. They found 
i = 44° and ; = 54° for AM 2322A and AM 2322B, respectively. 

The rotation curves and the spatial profiles of the Ha emis- 
sion and A 5735 continuum flux along the observed slit positions 
are presented in Figs.|4]and|5]for AM2322A and AM2322B, re- 
spectively. 

As can be seen in Fig. |4l AM 2322A does not have a well 
defined and symmetric rotation curve along to the observed slit po- 
sitions, indicating that the inclination of the galaxy must be quite 
low and that this object is actually being seen near face-on. Then, 
the real inclination of the galaxy must be lower than the one ob- 
tained from Ferreiro & Pastoriza (2004) and perhaps even from our 
own estimation. 

The heliocentric velocity of the main galaxy is taken to be the 
radial velocity measured at the nominal centre (continuum peak) 
along the slit position at PA=59° or v, = 3 739 km s"' . 

AM2322B shows a fairly symmetric rotation curve and for 
this galaxy we adopted a very simple approximation for the ob- 
served velocity distribution, assuming that the gas moves under a 
logarithmic gravitational potential, following circular orbits close 
to a plane P(i, t//(,), characterized by its inclination to the plane of 
the sky (/) and the position angle (PA) of the line of nodes ij/Q. This 
assumption results in an obs erved radial circular velocity v(r, ij/) in 
the plane of the sky given bv lBertola et alj ( Il99 ih : 

VqR cos(ifi - tf/o) sin(/) cos(0 



v(r,(/A)= V, + ■ 



^JR-r] + Rl cosHi) 



with 



[sin (i/r - ifio) + cos-(0 cos-((/^ - iflQ)], 



(1) 



(2) 



where Vs is the systemic velocity, R is the radius in the plane of the 
galaxy, and Vo and R^ are parameters that define the amplitude and 
shape of the curve. The fit of the rotation curve for this galaxy is 
shown in Fig.|5] 



4 Krabbe et al. 




Table 2. Kinematical parameters for AM Till^ 



Parameter 


AM2322B (PA=318°) 


in 


63 




318 


V, (km/s) 


3 376±6 


Vo (km/s) 


110±6 


Rc (kpc) 


-0.33±1391 



The above model for the rotation curve results in an heliocen- 
tric velocity of 3 376 km s"'. The observed radial velocities along 
the major axis are well represented by the model. The rotation curve 
is typical of spiral discs, rising shallowly and flattening at an ob- 
served amplitude of 110 km s"', with small (less than 25 kms"') 
deviations from the smooth rotational field, which are commonly 
obse rved in interacting galaxies. 

iDonzelli & Pastorizal ( 1 19971) estimated systemic velocities of 
3 680 and 3 424 kms"' for AM2322A and AM2322B, respec- 
tively, which agree within 2 % with our values. 

For AM 2322B, we can obtain an estimation of the dynamical 
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Figure 5. Observed Ho- and A 5735 flux (in units of lO" ergs cm"^ s"' ) and 
mean radial velocity as a function of apparent galactocentric distance along 
PA=318° for AM2322B. The velocity scale correspond to the observed 
and model values after subtraction of the systemic velocity of the galaxy, 
without correction by the inclination in the plane of the sky. 
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Table 3. Parameters used on the simulations 





AM 2322A 


AM 2322B 


Number of points in disk 


16384 


8192 


Disk mass 


0.6 


0.08 


Disk radial scale length 


0.7 


0.3 


Disk vertical scale thickness 


0.14 


0.0 


Reference radius K^ef 


1.7 


0.8 


Toomre Q at R^-^y 


1.5 


1.5 


Number of points in gas disk 


16384 


8192 


Gas disk mass 


0.06 


0.01 


Gas disk radial scale length 


1.0 


0.45 


Gas disk vertical scale thickness 


0.07 


0.01 


Toomre Q at Ki-ef 


1.2 


1.4 


Number of points in bulge 


1024 


1024 


Bulge mass 


0.5 


0.01 


Bulge radial scale length 


0.14 


0.05 


Number of points in spherical halo 


16384 


8192 


Halo mass 


2.4 


0.24 


Halo cutolf radius 


6.0 


2.0 


Halo core radius 


0.6 


0.4 



AM 2322B 

1 50 I ' ' ' 1 ' ' 




D [kpc] 



Figure 6. Circular' velocity curve of the initial model for AM 2322B (initial 
condition for the simulation). Rotation curves of the individual components 
are shown: spherical halo in orange, bulge in red, stellar disk in blue and 
gas disk in yellow. The total rotation curve is the continuous thick black 
line. Points with en'or bai's connected by thin black hne segments are the 
observed radial velocity values. 



Notes: Simulations were done in a system of units with G=l. Model units 
scales to physical ones such that: length unit is 3.5 kpc, unit velocity is 
262 km s"', unit mass is 5.586 X IO'^Mq and unit time is 13.062 Myr 

mass by assuming that the mass inside a certain radius is given by 
M{R) = RV-/G. Using the deprojected velocity amplitude of 110 
km s"' and a radius of 4 kpc its dynamical mass is 1.12 x 10'" A/©. 
It is important to emphasize that the maximum radius to which we 
can observe the gas in emission is quite certainly smaller than the 
total radius of the galaxies, so our estimate of the dynamical mass 
is a lower limit to the actual dynamical mass of the system. The 
estimation of the deprojected velocity is also dependent on the as- 
sumed inclination of the galaxies with respect to the plane of the 
sk y, so if we a ssume the inclination angle of ; = 54° derived 
by iFerreiro & P ast oriza (■2004.1 . the dynamical mass would be of 
1.36 X lO'^M©. 



4 NUMERICAL SIMULATIONS 

Aiming to reconstruct the history of the AM 2322-821 system and 
predict its evolution, we attempted to reproduce the interaction 
between AM2322A and AM2322B by running a series of N- 
body simulations. The simulations were carried out with the N- 
body/SPH code GADGET-2 developed by Springe! (2005j. Galax - 
ies were modeled following the prescription of IHemguisj (Il993h . 
where we included a gaseous disk component. 

The model parameters for AM 2322B were constrained from 
the observed moiphology and rotation curve presented in Sec. [3] 
The resulting simulated rotation curve is shown in Fig.|6] where we 
overplot the observed circular velocity data. 

In order to build a model for AM 2322A, we first attempted to 
measure its inclination through isophotal ellipse fitting (see Fig.O, 
but the results are inconclusive due to the small inclination angle, 
combined with tidal distortions experienced by the galaxy, which 
result in strong radial variations on the resulting isophotal posi- 
tion angle. The long-slit spectra crossing the disk through two dif- 
ferent position angles (see Figures [T] and |4), and none of those 
show clear signs of rotation, so the rotation curve cannot really 



be used to constrain the mass distribution of the galaxy as in the 
case of AM2322B. AM 2322 A is modeled as a face-on galaxy. 
A photometric mass estimate of 1.7 x 10'^ Mq^ is obtained us- 
ing a n absolute B magnitude Mg = -20.98 jFerreiro & Pastorizal 
2004) and assuming a mass-to-light ratio of ys = 4.4 M© Lg' from 
Faber & Gallaaher (1979) as valid for a Hubble type SA(r)c galaxy 
Ide Vaucouleurs et al.lll991t) . 

As usual in this kind of approach, to reproduce the dynamical 
and moiphological state of the AM 2322-821 system we have to 
solve a reverse problem of finding the orbit followed by the galax- 
ies from their observed properties, and this is not a fully determined 
problem, since the observational data do not provide all the neces- 
sary information. Therefore, in order to set up the initial conditions 
for the simulations, we calculate orbits that satisfy the requirements 
given by the observed radial velocity difference, testing different 
eccentricities, pericenter distances, and line of sight direction dis- 
tances. After that, based on the observed morphology and previous 
experience, we select a few orbits to simulate, from which the one 
that best fits the observed properties is selected. 

Different galaxy models were tested, and the parameters of the 
final ones are presented in Table |3] The AM 2322A model has a to- 
tal mass of 1.7 x 10"Mo, and the mass of the AM2322B model 
is 1.6 X 1O"'M0. The mass of each individual component (disk, 
bulge, halo, and gas) of both models is also given in the table. The 
disk inclination in the models were chosen to match the observa- 
tionally derived values (see Section |3] AM2322A is modeled as a 
face-on galaxy). The determination of AM2322B model inclina- 
tion required several test simulations, because the tidal distortions 
and warping induced by the interaction modifies the projected axial 
ratio: initial conditions with AM 2322B model inclined by / = 68° 
led to final morfologies that do not match with the observed one. 
The best results were obtained with AM 2322B initially inclined by 
i = 80°. A total of 75776 particles were used. 

After several runs, the orbit that best reproduces the obser- 
vational properties is a fast hyperbolic orbit, with an eccentricity 
e = 3.1 and perigalacticum of q = 10.5 kpc. The orbital plane is 
inclined to the plane of the sky by 37°, and intersects the later in a 
position angle of 80°. 
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Radial Velocity AM2322A, t = 91.27Myr 



Radial Velocity AM2322B, t=91 .27Myr 




• I f O « O 



• # i o 



Pos- [kpc] Pos. [kpc] 

Figure 7. Radial velocity curves of the models at the final stage, taken at the same position angles as the observed ones, for compaiison. Left and right panels 
are for AM 2322A and AM 2322B, respectively. Losanges are from the simulations, black dots with en'or bars are from observations. 



The dynamical mass of AM2322B, as obtained in Section|2] 
up to a radius of 4 kpc, amounted to 1.12 x IO'^Mq. At the final 
stage of the simulation, up to that radius, the model provided a mass 
of 1.1 X IO'^Mq (the initial model has a mass of 1.22 x W°Mq up 
to that radius). 

The simulation indicates that the spatial distance between both 
galaxies is 43.8 kpc (sky projected distance is 35 kpc). AM 2322A 
is closer to us than AM 2322B. 

Radial velocity curves of model galaxies AM2322A and 
AM2322B are presented in Figure |7] They are compared with the 
observed radial velocity data and show that the kinematics of the 
models is correct. 

Figure [8] shows the time evolution of the encounter. Time 
is shown in Myr in the upper right comer of each frame, with 
respect to the perigalacticum. Simulation starts ~ 120 Myr before 
perigalacticum. The situation that best reproduces the morphology 
and kinematics of the present stage of the AM 2322-821 system 
is at f = 91.27 Myr after perigalacticum. The general large scale 
morphology and kinematics agrees well with observations, within 
the resolution provided by the simulations. 



5 STELLAR POPULATION SYNTHESIS 

A detailed study of the star formation in minor merger galaxies is an 
important source of information not only on the age distribution of 
their stellar population components, but to better understand several 
aspects related to the interacting process, its elTect in the properties 
of the individual galaxies and their later evolution. The absorption 
features arising from the stellar component also affect to different 
degrees the measured intensity of the emission line in the spectrum 
of the gaseous component. This effect is more prominent in, but not 
restricted to, the Balmer lines, so the stellar population contribution 
must be subtracted from each spectra in order to study the physical 
properties of the gas. 

To investigate the star formation history of AM2322A 
and AM2322B , we use the stella r popu l ation s ynthesis code 
STARLI GHT llCid Fernandes et al.l 12004 l2005l : iMateus et all 
2006; A sari et al.1 l2()07h. This c ode is extensively discussed in 
Cid Fernandes et al. I( l2004l2005h . and is built upon computational 



techniques originally developed for empirical population synthe- 
sis with additional ingredients from evolutionary synthesis mod- 
els. The code fits an observed spectrum 0,i wit h a combination of 
Nj, sin gle stellar populations (SSPs) from the Bruzual & Chariot 
( I2OO3I) models. These models are based on a high-resolution library 
of observed stellar spectra, which allows for detailed spectral evolu- 
tion of the SSPs at a resolution of 3 A across the wavelength range 
of 3 200-9 500 A with a wide range of metallicities. We used th e 
Padova 1994 tracks as recommended bylBruzual & Chariot) (l2003h , 
with the initial mass function of Chabrier (Chabrier 2003) between 
0.1 and 100 Mq. Extinction is modeled by STA RLIGHT as due to 
foreground dust, using the reddening law of C ardelli et al. 
with Rv= 3.1 , and parametrized by the V-band extinction Ay. The 
SSPs used in this work cover fifteen ages, t = [0.001 , 0.003 , 0.005 , 
0.01 , 0.025 , 0.04 , 0. 1 , 0.3 , 0.6 , 0.9 , 1.4 , 2.5 , 5 , 1 1 , and 13] Gyr, 
and three metallicities, Z = [0.2 Zq, 1 Zq, and 2.5] Zq, summing 
up to 45 SSP components. Briefly, the code solves the follow ing 
equation for a model spectrum Mq dCid Fernandes et ^|2005|) : 



M, = M, 



A r.t 



G(v*,cr*) 



(3) 



were bj,^ is the reddened spectrum of the 7'th SSP normalized at 
Ao; = 10" "^'^•'"''•'o' is the reddening term; M,io is the synthetic flux 
at the normalisation wavelength; x is the population vector; igi de- 
notes the convolution operator; and G(v*, cr*) is the gaussian distri- 
bution used to model the line-of-sight stellar motions, it is centred 
at velocity with dispersion cr-,, . 

The fit is carried out with a simulated annealing plus Metropo- 
lis scheme, which search es for the minimum of the equation 
JCid Fernandes et al.ll2005h : 



(4) 



where emission lines and spurious features are masked out by fixing 
= 0. For more details on STARLIGHT see lCid Fernandes et al.l 
l l2005b . 

Prior to the modelling, the SSPs models were convolved to the 
same resolution of the observed spectra; the observed spectra were 
shifted to its rest-frame, corrected for foregroun d Galactic red den- 
ing of £(S - V) = 0.181 mag taken from .Schlegel et al.1 ( 1 19981) and 
normalized to /15870A. The error in 0,\ considered in the fitting 
was the continuum rms with a 5/A' > 10, where S /N is the signal- 
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Figure 8. Time evolution of the encounter. Time is sliown in Myr in tlie upper right comer of eacli frame, with respect to the orbital pericenter. The simulation 
begins at the upper left frame. Stars are shown in light colors, while the gas is shown in blue. Dark matter halo is not shown. Best fit to the current state of the 
AM 2322-821 system is shown in the last frame. 



to-noise ratio per A in the region around Aq = 5870 A. Also, the fit- 
ting was performed only in spectra with the presence of absorption 
lines. Measurement errors are still a problem in population synthe- 
sis. The most serious drawback of STARLIGHT, as it is, is that 
it does not provide error estimates on its para meters. The reliabil- 
ity of parame ter estimation was best studied in lCid Femandes et al.l 
( I2OO4I [200 5^ by means of simulations which fed the code with 
spectra generated with known parameters, add noise, and then ex- 
amined the coiTespondence between input and output values. They 
performed this kind of simulation for a. Nj,= 20 base and the main 
results of that study were: (1) in the absence of noise, the method 
recovers all components of x to a high degree of accuracy; (2) in the 



presence of noise, however, the individual output Xj fractions may 
deviate drastically from the input values. However, this problem 
can be circumvented by binning the stellar populations according to 
the fiux contributions. In addition to, the absorption features present 
in our spectra, i.e. Balmer lines, are mainly dependent on age rather 
than on metallicity. In Fig.[TT]is a histogram with the contribution 
of each stellar populations to the optical flux A 5 870 A of the spec- 
trum of the central region along the PA=59° , before the binning 
the stellar populations. About the same results were obtained for 
the other apertures. It can be noted in Fig. [TT]that the about the 
same contributions by age are obtained from different metallicities. 
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Figure 9. Stellar population synthesis for the central bin along the PA=59° and PA=318° slit positions, for AM 2322A (left) and AM 2322B (right), respec- 
tively. Top panel: spectmm coiTected for reddening and the synthesized spectrum (shifted by a constant). Bottom panel: pure emission spectrum. 
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Figure 10. Synthesis results in flux fractions as a function of the distance to the centre of the slit along PA=59° (left) and PA=28° (right) for AM 2322A. 



with a little increase of old population contribution when lower Z 
values are assumed (the metallicity-age degeneracy problem). 

Fig.|9]shows an example of the observed spectra corrected by 
reddening, the model stellar population spectra and the pure emis- 
sion spectra for AM 2322A and AM 2322B. The results of the syn- 
thesis are summarized in Tables |4] and |5] for the individual spatial 
bins in each galaxy, stated as the percentual contribution of each 
base element to the flux at i 5 870 A. 

Foflowing the prescription of Icid Femandes et al.l ( |2005|) . we 
defined a condensed population vector, by binning the stellar pop- 
ulations according to the flux contributions into very young, xyy (t 
< 1 X 10^ yr), young, xy (1 x 10^ < t ^ 5x 10* yr), intermediate-age, 
xi (5 X 10** < t < 2 X lO' yr) and old, xq (t > 2 x 10' yr) compo- 
nents. The same bins were used to represent the mass components 
of the population vector (myy, niy, mi and mo). The metallicity, an 
important parameter to characterize the stellar population content, 
is weighted by light fraction, and the results point out to a mean 
value around solar. The quality of the fitting result is measured by 



the parameters x'^ and adev. The latter gives the percentual mean 
deviation \0a - M,i\IOa over all fitted pixels, where 0,i and are 
the observed and model spectra, respectively. 

The spatial variation in the contribution of the stellar- 
population components are shown in [TO] and [T2] for the galaxies 
AM 2322A and AM 2322B, respectively. 

As can be seen in Fig.[TOl the stellar population of AM 2322A 
is heterogenous along both slit positions. It does not show age gra- 
dient: the very young xyy, intermediate xj and old populations xq 
contribute significantly to the optical flux at /1 5 870 A. In a mass- 
weighted context, the bulge is predominantly com posed of old pop- 
ulation, in agreement with the results obtained by ( IMac Arthur et al.l 
120091) . 

AM2322B is dominated, by the very young population com- 
ponent xvY. This component has a systematic variation along the 
slit, increasing outwards. Instead the fraction of old and interme- 
diate stellar population components is decreasing from the centre 
to the outer regions. The young star-formation episode ocurred in 
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Figure 11. Histogram for the central region along PA=59° showing the 
contiibution to the flux continuum (in percentage of each stellar age com- 
ponent for the 3 metallicity values considered). 
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Figure 12. Same as\W\ but along PA=318° of AM2322B 



this galaxy could be related with the perigalactic passage, that is, 
about 90 Myr after perigalacticum. The effects of the galaxy inter- 
action in the star formation history are much more conspicuous in 
the secondary galaxy than in the main galaxy. This is expected as 
the triggering of star formation events by tidal interactions appears 
to be more efficient in less massive systems. Similar results were 
also found for the pair AM 2306-721 (iKrabbe et al..200a) . 



6 GAS PHASE O/H ABUNDANCE 

Once the stellar population contribution has been determined, the 
underlying absoiption line spectrum can be subtracted to allow the 
measurement and analysis of the line emission from the gaseous 
component. The line intensities were measured using Gaussian line 
profile fitting on the pure emission spectra. We used the IRAF splot 



routine to fit the lines, with the associated eiTor being given as 
cr- = (tI^iji + crjj,,^, where iTcom and crii„g are the continuum rms 
and the Poisson error of the line flux, respectively. The relative er- 
ror in the flux estimates of the emission lines are lower than 20 
%. The residual extinction associated with the gaseous component 
for each spatial bin was calculated comparing t he observed Hy/H;3 
and H a /up ratios to the theoretical values in iHummer & StorevI 
( Il987h for an electron temperature of 10 000 K and a electron den- 
sity of 100 cm"\ The observed emission line intensities were then 
corrected by this residual extinction using the lHowarthl ( ll983h red- 
dening function. 

Determining accurate element abundances of the gas phase 
from optical spectra is critically dependent on measuring temper- 
ature sensitive line ratios, such as [Oiii](/i4959 + /1 5007)//! 4363. 
However, when doing spectroscopy of H ii regions with high metal- 
licity and/or low excitation, temperature sensitive lines such as 
[Oiii]/i4363 are found to be weak or unobservable, and empiri- 
cal indicators based on more easily measured line ratios have to be 
used to estimate metal abundances. 

Because we did not detect any temperature-sensitive emis- 
sion lines in our spectra, we used the i?23=([Oii]/i3727+[Oin] 
^495 9-F[Oiii]^5007)/H/3 vs. [O iii]^ 5007/[O ii]i 3727 i lMcGaughl 
Il99lh diagnostic diagram to estimate the metallicity of the gas 
in the star forming regions by comparing the observed values 
with a grid of photoionization models. The [N ii]A 6584/[o ii]A 3727 
line ratio was used to break th e R23 degen eracy for our 
data, as sugg ested by ( Ke wlev & Ellisorj I2OO8I) . We verify 
that all observed star forming regions of our galaxies have 
log([N ii]A 6584/[o ii]A 3727 > - 1 , thus they are placed in the upper 
R23 branch. 

The photoionization models were built using the code 
Cloudy/08.00 cFerland..200Z) . In each model, a stellar cluster was 
assumed as the ionizing source with the stellar ene rgy distribution 
obtained from the synthesis models Starburst99 jLeitherer et al.l 
[l999). We calculated models with metallicities of Z = 2.0,1.0, 
0.6, 0.4, and 0.2 Zq, ionization parameter log U = -1.0, -1.5, -2, 
-2.5, -3.0, and -3.5, with the stellar cluster having an upper stel- 
lar mass limit of M„p = 100 Mq, age of 2.5 Myr and formed by 
an instantaneous burs t. The solar va lue of 12 -1- log(0/H) = 8.69 
is taken from Alle nde Prieto et al. 

I (|200D. Similar models have 
been used to d escribe observatio nal data of H 11 regions in inter- 
acting (IKrabbe et alj|200l |2008|) as well as in i solated galaxies 
Dors & Copetti 2005h . The reader is referred to iDors & Copettil 
20061) or a full description of the models. Fig.[T3]shows the R23 vs. 



[O m]A 5007/[O ii]A 3727 diagram, with the observed values super- 
posed in the grid of computed models. Filled squares correspond to 
regions in AM2322A and open squares to regions in AM2322B. 
Figs. 1141 and 1151 show the estimated O/H abundance distribution 
as a function of the galactocentric distance for AM2322A and 
AM2322B, respectively. 

As we can see in Figs.ll4landll5l the oxygen abundances de- 
rived from the R23 vs. [O lll]/[0 11] diagram suggest an relatively 
homogenous O/H value of 8.6(±0.01) and 8.40(±0.05) across the 
disc of the galaxies AM 2322A and AM 2322B, respectively. 

The central O/H values of for the two galaxies are lower than 
those obtained for field galaxies o f the same luminosit y, in agree- 
ment with the results obtained by iKewlev et alj j2006! ) for galaxy 
pair members with small projected separations (s < 20 kpc). This 
result provides strong observational evidence for the hypothesis 
that the galaxy interactions create gas flows towards the central re- 
gions, caiiying less enriched gas from the outskirts of the galaxy 
into the central regions, mixing and homogenizing the chemical 
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Table 4. Stellar-population synthesis results for AM 2322A 
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(mag) 
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-7.98 


14.3 


22.8 


18.2 


40.6 


0.4 


2.3 


3.1 


94.2 


0.014 


1.4 


5.83 


0.15 


-7.56 


21.1 


26.3 


0.0 


49.9 


0.5 


1.8 


0.0 


97.7 


0.028 


0.1 


7.98 


0.00 


-7.14 


19.9 


13.7 


10.2 


56.6 


0.8 


1.8 


5.7 


91.6 


0.010 


0.1 


24.79 


0.00 


-6.72 


9.7 


25.9 


40.2 


21.2 


0.2 


4.0 


12.4 


83.4 


0.024 


0.1 


10.96 


0.00 


-5.88 


14.8 


0.0 


61.8 


20.7 


0.3 


0.0 


37.0 


62.7 


0.007 


0.1 


13.14 


0.42 


-5.46 


13.2 


25.6 


29.6 


29.4 


0.1 


3.8 


10.2 


85.9 


0.019 


0.0 
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0.48 
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37.4 


0.1 


3.3 


3.5 


93.1 


0.016 


0.1 


7.65 


0.24 


-4.62 


10.9 


38.9 


0.0 


46.8 


A 1 


z.z 


0.0 


97.7 


0.021 


0.1 


5.59 


0.25 


-4.20 


32.3 


33.7 


14.2 


16.8 


1.7 


10.7 


9.4 


78.2 


0.022 


0.1 


4.88 


0.34 


-3.78 


19.7 


25.5 


22.4 


31.0 


0.4 


2.0 


5.8 


91.8 


0.033 


0.1 


4.95 


0.12 


-3.36 


13.4 


6.9 


39.1 


39.5 


0.1 


0.7 


7.9 


91.2 


0.022 


0.1 


5.54 


0.18 


-2.52 


8.8 


8.5 


43.7 


37.6 


0.1 


0.9 


10.3 


88.8 


0.021 


1.2 


3.88 


0.38 


-2.10 


0.7 


0.0 


64.6 


33.7 


0.0 


0.0 


11.9 


88.1 


0.021 


0.1 


12.37 


0.05 


-1.68 


0.0 


17.8 


29.5 


51.8 


0.0 


0.8 


4.6 


94.5 


0.017 


0.2 


22.41 


0.21 


-1.26 


0.0 


28.4 


26.5 


43.7 


0.0 


1.5 


5.1 


93.4 


0.031 


0.1 


19.28 


0.51 


-0.84 


1.9 


25.8 


1.0 


69.7 


0.0 


1.3 


0.2 


98.5 


0.025 


0.0 


15.76 


0.78 


-0.42 


22.6 


28.8 


30.8 


17.2 


2.1 


5.4 


21.0 


71.6 


0.018 


0.1 


7.19 


1.40 
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0.1 
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Table 5. Stellar-population synthesis results for AM 2322B 
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composition of the interstellar medium dKewIev et alj|2006h . This 
is also si milar to the result p reviously found for the AM 2306- 
721 pair jKrabbe et al.ll200^ . where the companion galaxy also 
presented an oxygen abundance relatively homogeneous across the 
disc. 

We compare the metallicity gradients obtained for AM2322A 
and AM 2322B with the grad ients for the sample of close pairs de- 
rived by Kewlev et alj ^20 10) and for the isolated galaxies M 101, 
Milk Way, M83, and NGC300. The galactocentric distance is 
given in units of R/R25, where ^25 is the B-band isophote at 
a surface brightness of 25 mag arcsecond"^. The R25 adopted 
for AM2322A an d AM2322B are 13.5 kpc and 4.2 kpc, re- 
spectively. iKewlev et alJ l201Ci ) computed the gradients using the 
[Nn]^6584/[Oii]/i3727-0/H relation for the first three isolated 
galaxies and for their sample pair. The metallicity gradient in 
NGC 300 was computed using the emission line intensities of H 11 
regions obtained by Bresolin et al. (2009) and also using the method 
cited above. These gradients are shown in Fig. [16] For consis- 
tency, we recalculate the gradients in the AM2322A and AM2322B 
using the same method and found (— ^^^) = -0.13(±0.03) 

and i—^^—^) = -0.36(±0.03), respectively. As can be seen in 
Fig. [161 these gradients have about the same slope as close pairs 
(~ -0.25) and are flatter than the gradients in the isolated galax- 

ies ( 0.57) . These results are in agreement with those found by 

iKewlev et all (lolO) in galaxy pairs. We interpreted that this shal- 
lower gradients can be explained b y the action of inwa rd and out- 
ward radial flows of interstellar gas ( Krabbe et al .l2008h . In isolated 
spirals, this mechanism seems to be weak or non existent. 

N-body/SPH numerical simulations of equal-mass mergers 
predict that the radial metallicity gradients of the disk galaxies 
should flat ten shortly after th e first pericenter, due to radial mix- 
ing of gas teupke et alj2010l) . Although AM 2322-821 is a system 
composed by galaxies of quite different mass (A/prman /Mtfcourfan- ~ 
11 from the models in Section |4), our numerical simulations indi- 
cate that the current stage of the merger is be about 90 Myr after 
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Figure 13. The relation log(S23) vs. log([0 iii]/[0 11]) for the individual spa- 
tial bins in AM 2322B (filled squares) and AM 2322A (open tiiangles). The 
curves represent the photoionization models described in the text (dashed 
lines correspond to different values of the logarithm ionization parameter, 
solid lines to different gas metallicities). 



first passage, which agrees with the scenario proposed by the above 

authors. 

In our previous work dKrabbe et al.ll2008b . we found that for 
the AM 2306-721 pair only the companion galaxy presented an 
oxygen abundance relatively homogeneous across the disc, while 
the main galaxy showed a clear radial oxygen abundance gradient. 
Both systems, AM 2306-721 and AM 2322-821, are morphologi- 



12 Krabbe et al. 



2 4 6 8 

Distance (kpc) 

Figure 14. The 0/H abundance distribution as a function of the dis- 
tance to the center of the AM2322A estimated using the log(S23) vs. 
log([0 m]/[0 ii]). The solid line represents the linear fit of the data. 
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Figure 15. Same as Fig.[T4] but for AM2322B. 

cally similar. Each system is composed by a spiral and an irreg- 
ular galaxy. However, the mass ratios estimated from the numeri- 
cal simulations for AM 2306-721 is Mp,-j„,„rt/M secondary ~ 2 and for 
AM 2322-821 is about 11; and the numerical simulations predicted 
that the evolution of the encounter was more violent and slower 
in AM 2306-721. Then if we assumed that there was an oxygen 
abundance gradient before the encounter in the main galaxies in 
both systems, why the gradient was not destroyed or flattened by 
the gas inflow in the main galaxy in AM 2306-721? One impor- 
tant point is that star-formation rate in the nuclear region of the 
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Figure 16. Metallicity gradients for AM2322A (blue line) and AM2322B 
(orange line). For comparison, we show the metallicity gradients for the 
isolated galaxies M 101 (green line). Milk Way (red line), M 83 (filled black 
fine), NGC 300 (cyan line), and 8 interacting galaxies (dashed black lines). 

main galaxy of AM 2306-721 is about 10 times higher than the one 
derived for AM 2322-821. These highest star-formation rate could 
have increased the heavy element content in the central regions and 
thus could have produced steeper gradients, and in these case, high 
gas inflow rate would be necessary to destroy or flatt the abundance 
gradient in the main galaxy of AM 2306-721. 



7 CONCLUSIONS 

An observational study of the effects of the interaction in the 
kinematics, stellar population, and abundances of the galaxy pair 
AM 2322-821 was conducted. The data consist of long-slit spectra 
in the wavelength range of 3 350 to7 130 A obtained with the Gem- 
ini Multi-Object Spetrograph at Gemini South. The main findings 
are the following: 

(i) A fairly symmetric rotation curve for the companion 
(AM2322B) galaxy with a deprojected velocity amplitude of 110 
km s"' is obtained, and a dynamical mass of 1.1 - 1.3 x IO'^Mq 
within a radius of 4 kpc can be estimated using this deprojected 
velocity. 

(ii) Asymmetries in the radial velocity field were detected for 
companion, very likely due the interaction between the galaxies. 

(iii) In order to reconstruct the history of the AM 2322-821 sys- 
tem and to predict the evolution of the encounter, we modeled the 
interaction between AM2322A and AM2322B through numeri- 
cal N-body/hydrodynamical simulations. The orbit that best repro- 
duces the observational properties is found to be hyperbolic, with 
an eccentricity e = 3.1 and perigalacticum of q = 10.5 kpc; the 
current stage of the system would be about 90 Myr after peri- 
galacticum. 

(iv) The companion galaxy is dominated by a very young (t < 
1 X 10^ yr) population, with the fraction of this population to the 
total flux at A 5 870 A, increasing outwards in the galaxy disk. 
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(v) The stellar population of AM 2322A is heterogenous along 
of the slit positions observed. 

(vi) The oxygen abundance spatial profiles obtained for both 
galaxies are relatively homogeneous across the galaxy discs. The 
absence of an abundance gradient in these galaxies is interpreted as 
it having been destroyed by interaction-induced gas flows from the 
outer parts to the centre of the galaxy. 



ACKNOWLEDGMENTS 

Based on observations obtained at the Gemini Observatory, which 
is operated by the Association of Universities for Research in As- 
tronomy, Inc., under a cooperative agreement with the NSF on 
behalf of the Gemini partnership: the National Science Founda- 
tion (United States), the Science and Technology Facilities Council 
(United Kingdom), the National Research Council (Canada), CON- 
ICYT (Chile), the Australian Research Council (Australia), Min- 
isterio da Ciencia e Tecnologia (Brazil), and SECYT (Argentina). 

The authors would like to thank Volker Springel for provid- 
ing them with GADGET-2. I.Rogrigues also acknowledges Insti- 
tuto National de Pesquisas Espaciais - INPE/MCT, Brazil, for pro- 
viding computer time in one of its clusters to run the simulations 
presented here. A. C. Krabbe thanks the support of FAPESP, pro- 
cess 2010/01490-3.q 

We also thank Ms. Alene Alder-Rangel for editing the English 
in this manuscript. 



REFERENCES 

Allende Prieto, C, Lambert, D. L., & Asplund, M. 2001, ApJ, 
556, L63 

Arp, H. C, & Madore, B. F. 1987, A Catalogue of Southern Pe- 
culiar Galaxies and Associations 2 volume set, by Halton C. Arp 
and Barry F. Madore, pp. 504. ISBN 0521343364. Cambridge, 
UK: Cambridge University Press, May 1987 

Asari, N. v., Cid Femandes, R., Stasihska, G., Torres-Papaqui, 
J. P, Mateus, A., Sodre, L., Schoenell, W., & Gomes, J. M. 2007, 
MNRAS, 381, 263 

Barton, E. J., Geller, M. J., & Kenyon, S. J. 2000, ApJ, 530, 660 

Barton Gillespie, E., Geller, M. J., & Kenyon, S. J. 2003, ApJ, 
582, 668 

Bender, R., & Moellenhoff, C. 1987, A & A, 177, 71 

Bertola, R, Bettoni, D., Danziger, J., Sadler, E., Sparke, L., & de 

Zeeuw, T. 1991, ApJ, 373, 369 
Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 
Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 

245 

Chabrier, G. 2003, PASP 115, 763 

Cid Femandes, R., Gu, Q., Melnick, J., Terlevich, E., Terlevich, 
R., Kunth, D., Rodrigues Lacerda, R., & Joguet, B. 2004, MN- 
RAS, 355, 273 

Cid Femandes, R., Mateus, A., Sodre, L., Stasinska, G., & Gomes, 
J. M. 2005, MNRAS, 358, 363 

de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G., Jr., Buta, 
R. J., Paturel, G., & Fouque, R 1991, Third Reference Catalogue 
of Bright Galaxies. Volume III: Data for galaxies between 12* 
and 24*., by de Vaucouleurs, G.; de Vaucouleurs, A.; Corwin, 
H. G., Jr.; Buta, R. J.; Paturel, G.; Fouque, P.. Springer, New 
York, NY (USA), 1991. 

Donzelli, C. J., & Pastoriza, M. G. 1997, ApJS, III, 181 



Dors, O. L., Jr., & Copetti, M. V. R 2005, A&A, 437, 837 
Dors, O. L., Jr., & Copetti, M. V. R 2006, A&A, 452, 473 
Ellison, S. L., Patton, D. R., Simard, L., & McConnachie, A. W. 

2008, AJ, 135, 1877 
Faber, S. M., & Gallagher, J. S. 1979, ARA&A, 17, 135 
Ferland, G. J. 2002, Hazy, a brief introduction to Cloudy 96.03, 

Univ. Kentucky, Dept. Phys., Astron. internal report 
Ferreiro, D. L., & Pastoriza, M. G. 2004, A&A, 428, 837 
Ferreiro, D. L., Pastoriza, M. G., & Rickes, M. 2008, A & A , 481, 

645 

Friedh, D., Benz, W., & Kennicutt, R. 1994, ApJ, 430, L105 
Geller, M. J., Kenyon, S. J., Barton, E. J., Jarrett, T. H., & Kewley, 

L. J. 2006, AJ, 132, 2243 
Henry, R. B. C, & Worthey, G. 1999, PASP HI, 919 
Hemquist, L. 1993, ApJS, 86, 389 
Howarth, I. D. 1983, MNRAS, 203, 301 
Hummer, D. G., & Storey, P J. 1987, MNRAS, 224, 801 
Kennicutt, R. C, Jr., Roettiger, K. A., Keel, W. C, van der Hulst, 

J. M., & Hummel, E. 1987, AJ, 93, lOU 
Kewley, L. J., Dopita, M. A., Sutheriand, R. S., Heisler, C. A., & 

Trevena, J. 2001, ApJ, 556, 121 
Kewley, L. J., & Dopita, M. A. 2002, ApJS, 142, 35 
Kewley, L. J., Geller, M. J., & Barton, E. J. 2006, AJ, 131, 2004 
Kewley, L. J., & Elhson, S. L. 2008, ApJ, 681, 1183 
Kewley, L. J., Rupke, D., Jabran Hahid, H., Geller, M. J., & Bar- 
ton, E. J. 2010, ApJ,72I,L48 
Krabbe, A. C, Rembold, S. B., & Pastoriza, M. G. 2007, MN- 
RAS, 378, 569 

Krabbe, A. C, Pastoriza, M. G., Winge, C, Rodrigues, I., & Fer- 
reiro, D. L. 2008, MNRAS, 389, 1593 

Lambas, D. G., Tissera, P. B., Alonso, M. S., & Coldwell, G. 2003, 
MNRAS, 346, 1189 

Leitherer, C, et al. 1999, ApJS, 123, 3 

MacArthur, L. A., Gonzalez, J. J., & Courteau, S. 2009, MNRAS, 
395, 28 

Marquez, I., Masegosa, J., Moles, M., Varela, J.,Bettoni, D., & 

Galletta, G. 2003, Ap&SS, 284, 711 
Mateus, A., Sodre, L., Cid Fernandes, R., Stasihska, G., 

Schoenell, W., & Gomes, J. M. 2006, MNRAS, 370, 721 
McGaugh, S. S. I99I, ApJ, 380, 140 

Michel-Dansac, L., Lambas, D. G., Alonso, M. S., & Tissera, P. 

2008, MNRAS, 386, L82 
Nikolic, B., Cullen, H., & Alexander, P 2004, MNRAS, 355, 874 
Pagel, B. E. J., Simonson, E. A., Terlevich, R. J., & Edmunds, 

M. G. 1992, MNRAS, 255, 325 
Peeples, M. S., Pogge, R. W., & Stanek, K. Z. 2009, ApJ, 695, 

259 

Perez-Montero, E., & Contini, T. 2009, MNRAS, 398, 949 
Roy, J.-R., & Walsh, J. R. 1997, MNRAS, 288, 715 
Rupke, D. S. N., Veilleux, S., & Baker, A. J. 2008, ApJ, 674, 172 
Rupke, D. S. N., Kewley, L. J., & Barnes, J. E. 2010, ApJ, 710, 
L156 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 
525 

Sekiguchi, K., & Wolstencroft, R. D. 1992, MNRAS, 255, 581 
Vila-Costas, M. B., & Edmunds, M. G. 1992, MNRAS, 259, 121 
Sellwood, J. A., & Wilkinson, A. 1993, Reports on Progress in 
Physics, 56, 173 
Springel, V. 2005, MNRAS, 364, 1105 
Woods, D. R, & Geller, M. J. 2007, AJ, 134, 527 



100 



SS 75 - 



"5 50 



Very Young 



o 



25 - 




100 




llllll 




SS 75 - 



"5 50 
JO 

8 25 



100 

g 75 

o 



Young 



Intermediate 



o 
O 



50 - 



25 - 



100 



75 



Old 



■llllll I 



1 



